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Stability  Criteria  for  Elastic  Materials 

by 

K.N.  Sawyers  and  R.S.  Rivlin 

Center  for  the  Application  of  Mathematics 
Lehigh  University 
Bethlehem,  Pennsylvania 


Abstract 

Conditions  are  derived  which  are  necessary  for 
stability  of  inco npressible  elastic  materials.  These 
are  obtained  by  considering  small  shearing  displace¬ 
ments  superimposed  on  a  finitely  deformed  state  of 
the  material. 
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1.  Introduction. 

In  the  classical  elasticity  theory  for  infinitesimal 
deformations,  certain  restrictions  must  be  placed  on  the  elas¬ 
tic  moduli  to  ensure  that  the  material  be  stable.  For  example, 
in  the  case  when  the  material  is  isotropic  and  incompressible, 
the  shear  modulus  p  must  be  positive.  This  condition  en¬ 
sures  that  the  velocity  of  propagation  of  a  plane  sinusoidal 
shear  wave  be  real.  If  it  were  violated,  a  plane  spatially- 
sinusoidal,  shear  disturbance,  of  initially  infinitesimal 
amplitude,  would  increase  in  amplitude  with  time  even  though  no 
forces  are  applied  to  the  material.  Accordingly,  due  to  the 
presence  of  Brownian  fluctuations,  the  material  could  not  exist 
in  the  unde formed  state. 

Similarly,  in  the  case  of  incompressible  elastic  mater¬ 
ials  which  undergo  finite  deformations,  material  stability 
requires  that  the  strain-energy  function  be  such  that,  for 
each  state  of  pure  homogeneous  deformation  for  which  it  is 
assumed  to  be  valid,  the  velocities  of  plane  sinusoidal  shear 
waves  of  arbitrary  direction  of  propagation  and  infinitesimal 
amplitude  must  be  real. 

It  has  been  shown  that  this  implies  that  the  strain- 
energy  function  W  per  unit  volume  must  satisfy  the  conditions 
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where  A^  (A=l,2,3)  are  the  principal  extension  ratios  for 
the  pure  homogeneous  deformation,  and  W^g  (A,B  =  1,2) 

p 

denote  the  derivatives  9W/3I^  and  9  W/9I^3Ig  respectively 
of  W  with  respect  to  the  strain  invariants  I1  and  I2 
defined  by 


=  A?  +  A? 


I  =  A^X^+  A^A^ 

2  A2A3  A3A1  A1A2 


The  first  of  these  conditions,  which  was  obtained  by 
Ericksen  [1],  ensures  that  a  shear  ~ve,  polarized  in  the 
principal  direction  corresponding  to  the  extension  ratio  A^ 
and  propagated  in  an  arbitrary  direction  perpendicular  to 
this  principal  direction,  shall  have  a  real  velocity  of  propa 
gation.  The  second,  which  was  derived  by  Sawyers  and  Rivlin 
[2],  ensures  that  a  shear  wave,  polarized  in  the  plane  perpen 
dicular  to  the  principal  direction  corresponding  to  the  exten 
sion  ratio  A^  and  propagating  in  an  arbitrary  direction  in 
this  plane,  shall  have  a  real  velocity  of  propagation. 

In  the  present  paper,  it  is  shown  that  both  of  these  con 
ditions  can  be  derived  from  a  single  assumption.  We  consider 
that  an  infinitesimal  simple  shear  is  superposed  on  the  pure 
homogeneous  deformation.  The  plane  of  shear  is  a  principal 
plane  for  the  pure  homogeneous  deformation  and  the  direction 
of  shear  is  an  arbitrary  direction  in  this  plane.  The  assump 
tion  is  made  that  the  incremental  shear  modulus  is  positive. 


-  4  " 

2.  Basic  Equations 

We  consider  deformations  which  result  from  the  super¬ 
position  of  an  infinitesimal  deformation  on  a  finite  pure 
homogeneous  deformation.  Let  £  denote  the  coordinates  in 
a  rectangular  cartesian  system  x  of  a  generic  particle  of 
the  body  in  its  undeformed  state.  Let  x^  be  the  coordinates 
of  the  same  particle  in  the  deformed  state.  We  may  then  write 

xi  =  Xi  +  eui  C1) 

where  e  is  a  small  parameter  whose  squares  and  higher 
powers  may  be  neglected,  and 


xi  -  Mi 


*2  =  X252 


x3  -  x353 


The  A’s  in  Eqn.  2  are  the  extension  ratios  associated  with 
the  finite  deformation.  In  Eqn.  1,  the  displacements  u^ 
are,  in  general,  regarded  as  functions  of  X1 ,  X2>  X3  . 

The  finger  strain  tensor  .  ,  associated  with  the 

A 

resultant  deformation,  is  given  by 


3x. 


3Xk  ax, 

9?m 


Xi,k 


(3) 


*  The  usual  summation  convention  applies  to  repeated  sub¬ 
scripts.  The  comma  notation  ,i  is  used  to  denote  partial 
differentiation  with  respect  to  X^  . 
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We  write 


C- -  =  C. •  +  ec- - 
ij  iJ  iJ 


and  substitute  from  Eqns .  1  and  2  into  3  to  find 


C11  A1  *  C22  X2  *  C33  X3  * 


C..  =  0  for  i  f  j 


Al(<SilUj,l+  <SjlUi,l')  +  X2(<Si2Uj  ,2+  6j2Ui,2') 


+  A3(<Si3Uj,3+  6j3Ui,3} 


where  6^  is  the  Kronecker  delta 

Three  scaJar  invariants  T.  (A=l,2,5)  of  the  tensor 


are  given  by 


X1  Cii  ’  X2  2(CiiCjj  '  CijCij)  ’  1 3  det  ^ij5 


Writing 


IA  =  IA  +  eiA  ,  A  =  1,2,3 


and  substituting  from  Eqns.  4,  5  and  6  into  7,  we  find 


Xl+X2+X3  ’ 


2,2,  ,2,2,  -.2,2 


=  *2X3+X3X 


2  ■>  2  •>  2 


r+X^X*  ,  I-  =  \*K\ 

1  ±  cl  J  L  c. 


i  =  2  J  X^  u  +  X^u  +  X^u 

1  1  11,1  C.  d  9  £- 


i2  2  VX1^X2+X3^U1  ,1  +  X2^X3+X1^U2,2  +  X3  ^Xl+X2^  U3 , 3  1  ^1°') 


23  =  2X1X2X3^U1,1+  U2,2+  U3,3') 


From  Eqns.  8,  9  and  10,  it  follows  that  the  condition 
T  =  1  for  a  deformation  to  be  isochoric  yields 


X^XgX^  =  1  and  uk  k  »  0 


Let  and  FT  ,  respectively,  denote  the  components  of 

the  unit  vectors  normal  to  a  surface  in  the  material  before  and 
after  it  has  undergone  the  infinitesimal  superposed  deformation. 
Then  (see  [3],  for  example) 


N\  =  N-  +  en. 
ii  l 


where 


ni  NiNkN*uk,*  '  Nkuk,i 


From  Eqn.l  we  readily  obtain,  neglecting  terms  of 
higher  degree  in  s  than  the  first,  the  relation 


3_ 


-  eu 


k,j  3Xk 
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The  strain  energy  W  per  unit  volume  for  an  isotropic 
incompressible  elastic  material  depends  on  the  strain 
only  through  the  invariants  T"  and  ,  so  that 
W  =  W(T  T  ).  The  Cauchy  stress  tensor  IT. .  is  then  given  by 

1  9  d  I  J 


2[0V  -  l*2ClkCkj ) 


(15) 


where  P  is  an  arbitrary  hydrostatic  pressure  and  we  use 
the  notation 


WA  =  8W/3Ia  ,  A  =  1,2  (16) 

Noting  that  the  invariants  given  by  Eqn.  9  are  those  associa¬ 
ted  with  the  finite  deformation,  we  make  use  of  Eqn.  8  and 
expand  W  and  W2  in  a  power  series  about  I1  and  I2  . 
Thus, 

W  =  W  +  efW  i  +  W  i) 

1  1  1  11  1  12  2 J 

(17) 

W  =  W  +  £  fW  i  +  W  i  1 

2  2  1  12  1  22  2J 

where 


W 


A"  V1!,^ 


W 


32W 


AB 


3IA91B 


A,B  =  1,2  (18) 


I,, I, 


-  8  . 


Writing 


jij  *  zij  *  s°ij  •  p  ■  p  *  £p 


and  substituting  from  Eqns .  4,  8,  17  and  19  into  15  ,  we 


find 


Eij  *  2KV  IlWa3Cij-  W2CikCkj)  *  P{ij 


a.  .  -  2[(W,+  I,K„)c.  .  -  W-(C.,  c,  .+  c.,  C,  .) 
ij  1  ^  1  12  i]  2l  ik  kj  lk  kj' 


*  Kv  "ll^X  *  V12<-hh+  V  +  f2l) 


(W121l'*'  W2212:)CikCkjI  "  pSlj 


The  traction  acting  on  a  surface  whose  unit  normal 

is  ,  and  measured  per  unit  area  of  that  surface,  is  given 
by 


T.  =  E  •  .  N  . 
i  ij  3 


Writing 


T.  =  T.  +  et. 

li  l 


and  substituting  from  Eqns.  12  and  19  into  22  yields 


T.  =  E--N-  ,  t-  =  a..N.  +  E..n. 
l  in  1  i  i]  J  1 3  3 


9 


In  the  absence  of  body  forces  the  equation  of 


equilibrium  is 


=  0 


Substitution  from  Eqns.  14  and  19  into  this  yields  the 
two  equations 


Sij,j  "  0  and  °i j  >  j  '  Uk,jZij,k  0  (26) 


We  now  assume  that  the  superposed  infinitesimal  deforma¬ 
tion  is  a  simple  shear  in  the  plane  perpendicular  to  the  3- 
direction. 

The  displacement  u^  associated  with  this  deformation  is 
given  by 


u^  =  k(X2cos0  -  X1sin6)cos6 

u2  =  k(X2cos0  -  X-sip0)sin0  (27) 


where  k  is  the  amount  of  shear  and  0  is  the  angle  between 
the  direction  of  shear  and  the  x-^-axis,  as  shown  in  Fig.l. 

We  remark  that  the  displacement  described  in  Eqn.  27  satisfies 
the  incompressibility  requirement  given  by  Eqn. 11. 

If  k  is  positive,  then  the  direction  of  shear  is  in 
the  direction  of  the  unit  vector  K  with  components 


-  cose  ,  K2  =  sine  ,  =  0  (28) 

Let  N  be  a  unit  vector  in  the  plane  of  shear  perpendicular 
to  K  ,  the  components  of  N  being  given  by 

N1  =  -  sine  ,  N2  =  cos0  ,  =  0  (29) 


We  consider  a  material  plane  in  the  body  which  has  undergone 


the  pure  homogeneous  deformation  described  by  Eqn.  2.  The 
unit  normal  N  to  this  plane  after  the  superposition  of  the 
infinitesimal  shear  is  given  by  Eqn.  12,  where,  from 
Eqns.  13,  27  and  29, 


n  =  0 


c 


We  substitute  from  Eqns.  27  into  6  and  10  and  obtain 

2  2 
n  =  -2AnK  sin6cos8  ,  c^2  =  2A2k  sihdcosQ  , 

C12  =  K C^2COS^6  ‘  Aisin20^  »  c23  =  C31  =  C33  =  0  *  (- 

i2  =  X^  =  2 X^k  sinOcosO (x|- X^) 

We  substitute  from  Eqns.  5  and  9  in  Eqn.  20  and  obtain 
the  components  of  the  Cauchy  stress  associated  with  the 
finite  homogeneous  deformation. 

in  -  ni  h ♦  '4 + 

J22  ■  2A2  {*1  +  <*3  + 

-  2X2,  {w,  ♦  (X?  * 


Similarly,  substitution  from  Eqns.  5,  9  and  31  in  Eqn.  21 
gives  the  additional  stresses  that  arise  from  the  deforma¬ 
tion  described  by  Eqn.  27.  Thus, 


4«A^sin0cos0  +  (A^-A^jw 


11+2X23W12 


+  X3W23  *  X2(W12*X3W22>}] 


4>cA^sinecose[wi+  A^-C^-A2)  / 


\l+W3W12 


.  x‘w22  ♦  a®(»12+  a23w22 


>}]■ 


o33  -  4cA3sin9cos0(A|-A3)[w2+W11.(A3+A2+A|)W] 

*•  X3W22(Ai+A2)]  '  P 

°12  =  2K{A2cos2e-A3sin2e)  Wj^+A^J 


°23  =  °31  =  0 


Bearing  in  mind  that  A^Ag.A  ,  <  and  0  are  constants, 

it  follows  that  the  stresses  given  by  Eqns.  32  and  33  satisfy 

the  equilibrium  conditions  of  Eqn.  26  provided  P  and  p  are 

constants . 

The  traction  acting  on  a  plane  with  unit  normal  N  can 
be  obtained  by  substituting  from  Eqn.  29  into  24.  By  making 
use  of  Eqn.  30  and  the  fact  that  certain  components  of  the 
stresses  vanish  (see  Eqns.  32  and  33),  we  obtain 

Tx  =  -  sine  Zi;L  ,  T2  =  cos 0  E22  ,  T3  =  0  (34) 


and 


\x;;  jcanssvsrr  ~rxc 


-  13 


=  -sinOa.^  +  cos0a12 


t2  =  -sin6a12  +  cos0a22 


t3  -  0 


Let  S’  denote  the  component  of  the  traction  vector  in 


the  shear  direction  K  .  Then,  from  Eqn.  23, 


S  =  T*  K  =  (T .  +  £ t . ) K - 

v  i  y  3 


Writing 


S  =  S  +  es 


and  substituting  from  Eqns .  28,  34,  and  35  into  36,  we  obtain 


S  =  sin0cose  (Egg-E-jj  ) 


2  2 

s  =  sin0cos0  (a22"°n  x)  +  ai2^cos  e"s^-n  0) 


Expressions  for  S  and  s  in  terms  of  material  properties  can 


be  obtained  by  substituting  from  Eqns.  32  and  33  into  38. 


Introducing  the  notation 


2(W1+Apf2) 


mi 


4.  Criteria  for  Material  Stability 

In  order  that  the  material  be  stable  for  given  values  of 
Xi , Xg >  and  X^  ,  it  is  necessary  that  the  superimposed  -ac¬ 
tion  have  a  positive  component  in  the  direction  of  the  super¬ 
imposed  shear  deformation,  i.e.,  the  shear  modulus  must  be 
positive . 

According  to  Eqns.  36  and  37,  s  is  the  extra  shear 
traction  that  is  required  to  produce  the  deformation  given  by 
Eqn.  27.  Since  k  is  the  amount  of  shear,  we  define  the  shear 
modulus  y  by 

V  =  ~/<  (41) 

Then  from  Eqn.  4G  we  obtain 

y  =  M(A2sin29+A2cos20) 

+  (X1-X2)2m  sin20cos20  (42) 

and  it  follows  that  a  necessary  condition  for  material  stab¬ 
ility  is  y  >  0  for  all  possible  values  of  0  .  In  particu- 

IT  3tt 

lar,  if  0  assumes  any  of  the  values  (0  ^ir,-^-)  then  y  >  0 
if  and  only  if 

M  >  0  (43) 

For  the  following  discussion  we  assume  that  Eqn.  43  is 
satisfied. 


Then  if  A1  =  ,  it  follows  that  y  >  0  regardless 

of  the  value  of  m  .  But  if  A  1  Ag  and  if 

m(A1-A2)2  =  -  MCA1+A2)2  (44) 

then  Eqn.  42  becomes 

y  =  M(A1sin2e-A2ccs20)2  (45) 

and  it  follows  that  y  is  non-negative  for  all  values  of 
0  and  vanishes  if 


tan  0  =  ±  /A27'A1 


(46) 


It  is  clear  from  Eqn.  42  that  for  any  fixed  value  of  0  , 

(other  than  those  for  which  either  sin  0  or  cos  0  vanish) 

y  increases  as  m(A1~A2)  increases  and  decreases  as 
2 

m(Ai~A2)  decreases.  Thus,  the  condition  y  >  0  for  all 
values  of  0  is  satisfied  if 


m  > 


M 


<yy2 

y-y2 


(47) 


However,  if  m(A1-A2)2  <  -  M(A  +A2)2  ,  then  y  <  0  for  0 
given  by  Eqn.  46,  and  we  conclude  that  the  material  is  in¬ 
herently  unstable. 

We  have  derived  two  conditions  that  are  necessary  for 
material  stability.  By  substitution  from  Eqn.  39  in  43  and 
47,  these  are  found  to  be 


-..^SRg&tsraBII&w”' ' :  —  ■SB' 
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W1  +  X3W2  =■  0 


W11  *  2X3W12  +  X3W22  > 


1  VX3W2 

2  (x,-x2)2 


Further  conditions  analogous  to  those  in  Eqn.  48  can  be 
derived  by  considering  superposed  shearing  deformations  similar 
to  that  of  Eqn.  27,  for  which  the  shear  planes  are  the  23  and 
31  planes,  respectively.  An  analysis  which  parallels  that 
given  above  then  yields  necessary  conditions  for  material  stab¬ 
ility  in  the  forms 


W1  +  X1W2  >  ° 


,  t  ,  vx2!*, 
W11  +  2AIW12  +  XIW22  >  -  I  ^7f2 


W1  *  *iW2  >  0 


2  1,  1  VX2W2 

W11  *  2X2W12  +  X2W22  »  -  I 


The  conditions  given  in  Eqns .  48,  49  and  50  may  be  re¬ 
cast  in  the  succinct  forms  shown  in  the  Introduction  by 
employing  Eqns.  9  and  11. 
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Figure  Caption 


Figure  1.  Geometry  of  a  small  superimposed  shearing  deforma¬ 


tion. 


